Abstract 



We present an analytical description of the motion in the singular loga- 
rithmic potential of the form = In \J x\/b 2 + x|, a potential which plays 
an important role in the modeling of triaxial systems (like elliptical galax- 
ies) or bars in the centers of galaxy disks. In order to obtain information 
about the motion near the singularity, we resort to McGehee -type trans- 
formations and regularize the vector field. In the axis-symmetric case 
(b = 1), we offer a complete description the global dynamics. In the non 
axis-symmetric case (6 < 1), we prove that all orbits, with the exception of 
a negligible set, are centrophobic and retrieve numerically partial aspects 
of the orbital structure. 

Key words: non axis-symmetric logarithmic potential, regularization, 
dynamics near the singularity. 
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1 Introduction 

The non axis-symmetric logarithmic potential plays an important role in galaxy 
dynamics. In the three-dimensional space, the potential: 



(where (xx,X2) are the usual cylindrical coordinates x\ = z, x 2 = R) models an 
elliptical galaxy with a dense core of radius R c and with the additional property 
of having a flat rotation curve at large radii In the 2D space, the logarithmic 
potential can describe other non axis-symmetric components of galaxies, such as 
bars in the centers of galaxy disks. 

The study of the orbital structure of the logarithmic potential was initially 
motivated by the need to construct self-consistent models of galaxies [HI . Nu- 
merical experiments have proven to be very useful in revealing the rich orbital 
structure of this potential, including the major orbit families, the resonances and 
the stochastic orbits. For example, in the axis-symmetric case (b = 1), the log- 
arithmic potential has been shown to admit only loop orbits, which are regular 
and avoid the origin. The non axis-symmetric {b < 1) potential admits two major 
families of orbits: box (for p = V R 2 + z 2 <C R c ) and loop orbits (for p ^> R c ) 
(DP' CI])- Note that this behavior can be also be retrieved from the two simple 
analytical approximations of the potential, one as a sum of two oscillators very 
close to the origin and the other one as ~ ln(p) at large distances 0. 

An interesting change in the orbital behavior has been discovered when the po- 
tential becomes singular. In their numerical study, Miralda-Escude & Schwarzschild 
[TO] have found that as R c — > 0, a larger fraction of regular box orbits becomes 
irregular or "box-like" (i.e. they will admit fewer integrals of motion than the 
number of spatial dimensions), with the end result that, in the limit R c = 0, 
all box orbits are irregular. The general interpretation of this result is that the 
scattering by the singularity renders the box orbits unstable, a similar behavior 



<S> = hn(R c 2 + 




(1.1) 
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to that observed in systems which contain a central black hole jSj. In addition, 
the singular logarithmic potential admits several families of minor orbits, i.e., 
resonances in terms of the X2-X\ frequency ratio: the banana (2:1), fish (3:2) and 
pretzel (4:3) orbits jTOj. Numerical studies have also revealed the existence of 
some stochastic orbits in a narrow region near the singularity [TK]). This 

has led some to suggest a link between the scattering by the singularity and 
the transition to chaos [3] , although no rigorous proof has been given to date in 
support of this hypothesis. 

The present investigation is an analytical approach to the study of the orbital 
dynamics (including the behavior near the singularity) in the case of the singular 
logarithmic potential: 



in both the axis-symmetric and non axis-symmetric cases. The aim of our work is 
not only to complement the previous numerical studies performed on this subject, 
but also to offer a theoretical basis for interpreting their results. 

The main difficulty in investigating the system analytically is due to the pres- 
ence of the singularity in the origin, which creates a discontinuity in the equations 
of motion. This problem motivates the introduction of a change of coordinates 
that regularizes the equations of motion. For this, we resort to McGehee-type 
transformations |2j , a technique that is frequently used in celestial mechanics for 
the study of singularities in the n-body problem (n = 1, 2, 3, ...). The underlying 
idea behind this technique is to transform the equations of motions and the time, 
such that the singularity is "blown- up" into a non-trivial manifold (in our case 
a torus). By studying the characteristics of the flow on this manifold, one can 
extrapolate the information (by continuity with respect to the initial data) about 
the orbital dynamics around the singularity |13j . 

The paper is organized as follows: Sections 2 and 3 provide a brief description 
of our system, in terms of the equations of motion and the conservation of energy. 
In Section 4, we remove the singularity by regularizing the equations of motion. 
In Section 5, we provide a description of the collision manifold and the zero 
velocity manifold, two abstract surfaces on which we can visualize the properties 
of the flow close to the singularity and at the maximum distance from the source 
allowed for a given energy, respectively. Section 6 presents the complete global 
dynamics in the axis-symmetric case (b = 1). In Section 7, we extend the analysis 
for the case of the non axis-symmetric potential (b ^ 1), in the restricted limit 
in which the anysotropy is small. In this case, we prove theoretically that the 
majority of orbits in this case are centrophobic (that is, they avoid the origin) — 
a result that has been originally discovered in the numerical study of Miralda- 
Escude Sz Schwarzschild jTUj — and we discuss the orbital structure in terms of 
orbits which preserve or change the sign of their angular momentum. Finally, in 
Section 8, we show how several families of resonances can be retrieved through 
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the numerical integration of the new system of equations. 



2 The Equations of Motion 

The non axis-symmetric logarithmic problem is a one-parameter Hamiltonian 
system with two degrees of freedom. The anisotropic logarithmic singular poten- 
tial (|l-2j) determines a conservative system with a preserved Hamiltonian given 
by: 

H(x u x 2 , yi ,y 2 ) = ^( Vl 2 + y 2 2 ) + H^jjr + (2.1) 

where (21,22) £ — {(0;0)} are the generalized coordinates and (2/1,2/2) G R 
are the momenta. 

We start by writing the equations of motion in a form that contains the 
anisotropy in the kinetic term rather than in the potential. Thus, we substitute 
1i — fr, <?2 = 2 2 , and Pi = yi/b 2 , p 2 = y 2 . Introducing the (standard) notations 
for the generalized vector coordinate q = (qi,q 2 ), the generalized momenta p = 

(pi,p 2 ), the anisotropy parameter a = — and the mass matrix 



M 

the Hamiltonian writes: 



[L 

1 



^(q,P) = ^P T Mp + ln|q| (2.2) 

and the anisotropic logarithmic problem is given as the first order system of 
ordinary differential equations: 

= Mp 

(2-3) 

w 

When jj, = 1, the above system describes the motion in the axis-symmetric 
logarithmic potential. From the point of view of Hamiltonian mechanics this 
case is completely integrable, as we have the two integrals of motion given by the 
conservation of the total energy and the total angular momentum. However, we 
point out that for |q| — > the dynamical behavior becomes unknown since the 
vector field (q, p) ceases to exist. 

Let K be the kinetic energy, 

K = ^p r Mp, (2.4) 

and V the potential, 

V = \n\q\. (2.5) 
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The total energy E is then: 



H(q,p) = K + V. (2.6) 

Since the system (J2.3j) is Hamiltonian the total energy is conserved. That is 
H(q, p) is constant along the solution curves of ()2.3j) and consequently the level 
sets of H (q, p) are invariant under the flow. If h is a certain energy constant (i.e. 
H((l,p) = h =constant) then the level set H(q, p) _1 (/z) is a three-dimensional 
surface usually called an energy surface, which we denote T, h . 

3 Topological Description of the Energy Sur- 
faces 

We want to consider in more detail the topology of various surfaces S^, including 
the orbits near the singularity. Let us fix h G R. Using a technique similar to 
McGehee's jU] we introduce the change of variables: 



re r 2 s 
iu 

r 



(3.1) 



where r > 0, s is a point on the unit circle S 1 and u G M 2 . Our transformation 



is a diffeomorphism from R 2 \ {(0, 0)} x R 2 to (0, oo) x S 1 



(where S 1 is the 



unit circle) and can be understood as a passing to some unorthodox kind of polar 
coordinates. 

The system f!2.3|) becomes: 



r 2 + 2 
1 i 
s = — 



e?s T Iu 



u 



Mu 
1 



r + 2 



- (s T Mu) ■ s 
s T Mu) ■ u - s 



(3.2) 



and the conservation of energy transforms to: 

-u T Mu + r 2 In r - 1 = h ■ r 2 . 



(3.3) 



Notice that the new system (j3.2j) is analytic on the open manifold (0, oo) x 
S 1 x M. 2 and that the regions of motion are constrained by the energy relation 
(13. 3 j) . More precisely, since the kinetic term in ()3.3|) is positive, we have that for 
a fixed level of energy h: 



hr 2 — r 2 lnr + 1 > 0. 



(3.4) 
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Solving the above relation, it follows that < r < R ma x with R m ax = Rmaxih). 
In other words the motion is always bounded for any fixed level of energy h. 
Also, we can say that the energy surface projects onto a disc of radius Rmax 
in R 2 \ (0,0). Since 



it follows that along the r = R m ax boundary we have u = 0, that is the kinetic 
term cancels. For this reason the curve 



is called the oval of zero velocity in S/j. In physical space, this curve represents 
the outmost boundary of the orbital structure allowed for a given energy. We 
will denote this curve by Z. 

PROPOSITION 3.1 For any fixed level of energy h the energy surface S/j is 
diffeomorphic to an open solid torus (i.e. a solid torus minus its boundary). 

Proof The proof follows closely the proof of Proposition 1.1 in 0j. Let B be the 
interior of the ellipse 



in the plane. Clearly B is topologically equivalent to S 1 x (0, Rmax] (r = 
corresponds to the missing boundary of the ellipse). Then we define F : — ► 
S 1 x B by 

F(r,s,u) = (s,u). 
F is just the required diffeomorphism. 

4 Regular izat ion of the Vector Field 

The objective of this section is twofold: to extend the energy relation to r = 0, and 
to regularize system (|3.3|) such that the new vector field becomes differentiable 
over the entire interval r G [0, R max ]- 

The energy relation ()3.3|) is ill defined at the singularity r = 0. However, the 
function r 2 In r can be continuously extended for r = by: 



-u T Mu = hr 2 — r 2 In r + 1 



(3.5) 




(3.6) 



-u r Mu = 1 

2 




(4.1) 
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The extended function (|4.1|) is also differentiable over its domain. Therefore we 
are able to extend the energy manifold by 




(4.2) 



for all r E [0,R max \. 

Now we will introduce a sequence of transformations of the system (J3.2J) such 
that the new system will have no singularity at r = 0. Following the technique 
introduced by McGehee pj, we want to paste an invariant manifold onto the 
phase space, such that we close the open solid torus that bounds the motion, by 
including its boundary r = 0. Also, in order to preserve the continuity of the flow 
with respect to the initial data, we have to ensure that the transformed system 
has a differentiable vector field. 

1 „ 

We implement a change of the time variable through da = — e^r dt. This will 
have the effect of decreasing the rate of the time intervals near the singularity. 
Expressed in the new time derivative the system (13.2)1 becomes: 



We note that by the above sequence of reparametrizations we have obtained an 
analytic vector field for (r, s, u) G [0, Rmax] x S 1 x M 2 , which is coupled with a 
differentiable integral relation given by ()4.2j) . 

It follows from (|4.2j) that each energy surface meets the boundary r = 
along a submanifold given by 



This manifold, let us call it A, is diffeomorphic with a two-dimensional torus, 
which we shall call the collision manifold. This fictitious torus has no meaning 
in the physical space, since the motion ceases to exist in the origin. However, 
the behavior of the orbits near the singularity can be extrapolated from similar 
properties of the flow on the collision manifold. In this sense, we regard the 
differential system of equations as a vector field on the manifold, its solutions 
representing the flow. 

Note that A is independent of the total energy and therefore we can say that 
the time transformations we have applied have the effect of pasting an invariant 
boundary onto each E^. Over this boundary the vector field is given by 




(4.3) 




(4.4) 



s = (s T Mu) ■ s - Mu 
u = 



(4.5) 
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The solid and now compact torus E^ bounds all the orbits and by investigating 
its properties, one can obtain a global picture of the motion, including the motion 
near the singularity. We emphasize that the orbits of system (|2.3J) are the same 
as the orbits of system (JH3J), only the parametrization is different. Therefore, 
any results concerning the solutions of the first system can be seen as results for 
the solutions of the second, as long as one is aware of the fact that the rate at 
which solutions move along the orbits is different. 

Using the energy integral we can further reduce the dimension of the system. 
For this, we express the coordinates (s, u) in terms of the new angle coordinates 
9 and ip\ 

s = (cos 9, sin 9) 

u = \/2(hr 2 — f(r) + 1)( cost/', svaip) (^-6) 

Also, in order to further simplify the system, we perform a similar time parametriza- 
tion which incorporates a part of the radial dependence in a new time variable: 

ds = ^/2(hr 2 - f(r) + l)da. (4.7) 

In these variables, the system (|4.Hjl becomes a first order system for (r, 9, ip) G 
[0, Rmax] x S 1 x S 1 with a differentiable vector field: 

2 — ^ j2(hr 2 — f(r) + l)(y / 7tcos6>cos'?/> + sin 9 sinip) 

9 = 2(hr 2 - f(r) + 1) (^/Ji cos ip sin 9 - sin^cosfl) ( 4 - 8 ) 
= r 2 (cos -ip sin 9 — y/Jl sin ip cos 9) . 

Before discussing different aspects of the dynamics of the above system, we 
make two more observations: 

OBSERVATION (1) Up to this point, no reference was made to the angular 
momentum integral. It is well known that in the axis-symmetric case (/i = 1), the 
singular logarithmic system admits, besides the total energy, a second conserved 
quantity, namely the angular momentum C := q2Pi — ?iP2- in terms of (r,9,ifj) 
this reads: 

C(s) = e- 1/r2 2[hr 2 - f(r) + 1] [^sintfcos^ - cos 9 sin ip], (4.9) 

where s is the final time variable introduced in (|4.7)l . 

Like the energy relation, C(s) is ill defined at the singularity. In an analogous 
way, we extend the function e _1//r continuously at r = by defining the function: 

,-i/r 2 if r >0 
9(r) = { n .. . (4.10) 

it r = 
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The angular momentum is now well defined and differential) le for all r > 0. For 
the axis-symmetric case (/i = 1), the law of conservation of angular momentum 
ensures that: 

as 

and therefore we have the integral relation: 

C(s) = g{r)2[hr 2 - f(r) + 1] sm(9 - V) = const. (4.12) 

In the anisotropic case (/j / 1) the above symmetry is lost. The anisotropy 
is responsible for the much more complicated dynamics and, eventually, for the 
existence of the chaotic motion. For later purposes, we write the variation of 
angular momentum 

dC 

— = (l-^p lP2 , (4.13) 
in terms of the new variables and time: 

^-(e-y^e) = ^—^e- 1/r2 2[hr 2 - f{r) + 1] 3/2 sin $ cos ij). (4.14) 



OBSERVATION (2) The orbital dynamics of the system (J4.8)) does not depend 
on the level of energy h. More precisely, each different but fixed h gives rise to 
the same qualitative phase portrait. The only change is in the value of R max , as 
Rmax depends directly on h. Therefore, without losing generality, and in order to 
simplify the calculations, we choose from now on to work with h = 0. 



5 The Collision Manifold and the Zero Velocity 
Manifold 

The flow on the invariant collision manifold A is given by imposing the restriction 
r = to the system ()4.8j) : 

9 = 2(J~jl cos ip sin 6 — sin ip cos 6) 

(5.1) 

7p = 0. ^ J 

We obtain a family of solutions ip — ipQ — const, whereas the vector field van- 
ishes along the deformed circles {(9, ip) \ i> = ^/Jicos1p sm9 — smipQCOs9 = 0} 
Notice that for u — 1, the curves of equilibria transform into circles given by 
{(9, ip) | ip — ip , 9 = ip , 9 = n + ip }. For the general case \i ^ 1, let us denote 
by C + the equilibrium curve that passes through 9 = ip = and by C~ the curve 
that passes through 9 = ip = n. It is immediate that C + is a repeller and C~ is 
an attractor. The dynamical behavior on the collision manifold A is illustrated 
in Figure 1. Since ip = 0, the flow follows the parallel lines ip = ipo = const. 
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Similarly, for r = R ma x we obtain another invariant manifold = for 
(JHHJ), namely the manifold corresponding to the oval of zero velocity. Recall that 
Rmax is the value which cancels hr 2 + 1 — /(r) = 1 — f(r), (h is set to be zero). 
The flow on Q is given by: 



The dynamical behavior on the zero-velocity manifold is similar to that on the 
collision manifold: the torus (r, 9, ip) G R ma x x S 1 x S 1 that represents Qo is 
covered by orbits parallel to 9 = 9q = const. There are again, two skewed circles 
of equilibria, {(9,ip) \ 9 = 9 , cos^sin^ — ^/Ji sin ip cos 9 = 0, }. Denoting by V + 
the curve that passes through (0, 0) and by V~ the curve that passes through 
(0, 7i"), it follows that V + is an attractor and V~ is a repeller (see Figure 2). 

6 Global Dynamics of the Axis- Symmetric Sys- 
tem (/i = 1) 

We return now the full system ()4.8|) and, taking into account that (at least in a 
prime analysis) the anisotropy is given by values of \x close to 1 but greater than 
1, we define the parameter e := ^fji — 1 > 0, which we will later treat as a small 
perturbation to the isotropic system. We also choose to work with the relative 
angle (p :— 9 — ip, instead of the angle 9. Then, written in terms of e and (r, tp, ip) 
and neglecting the terms of order e 2 and higher, the system ()4.8|) becomes: 



In the absence of the small perturbation (e = 0), that is in the axis-symmetric 
case 1^ — 1, the above system reduces to: 



with (r,<p,i/;) on the solid torus [0, Rmax] x S 1 x S 1 . The equilibrium solutions 
form two circles along ip = and <p = it and are given by (0, 0, V'o) an d (0, 7r, ip Q ) 
where ipo can be any value in [0, 2tt]. 




(5.2) 




(6.1) 




(6.2) 
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There are four invariant manifolds: 

- the collision manifold (at r = 0), on which the dynamics is given by: 




(6.3) 



- the zero velocity manifold (at r = Rmax), on which we have: 

= -RL» sin ¥> 
dip 

or just — = -1; 

dip 

- the "simp = 0" manifolds (when ip = or ip = n), with: 



(6.4) 



, _ 2 -2(l-/(r)),if^ = 
2(l-/(r)), if¥> = 7r. 



(6.5) 



r 2 + 2 



Recall that in the unperturbed case the angular momentum is conserved. In 
this case (see Observation 1): 

C(s) = g(r)2(l - f(r)) simp = C = const, (6.6) 

for all r e [0, R max \- It is easy to see that 

C = •<=>- [r = or r = -R ma x or sin cp — 0] . (6.7) 

In other words, the angular momentum is null (i.e. the motion is rectilinear) 
if and only if the orbits are either on the collision manifold, or on the zero 
velocity manifold, or connecting the two. From this perspective, we can denote 
the surfaces "simp = 0" as manifolds of zero angular momentum. Thus we have: 

PROPOSITION 6.1 In the axis-symmetric singular logarithmic problem the only 
orbits reaching the collision are the rectilinear ones. 

PROPOSITION 6.2 In the axis-symmetric singular logarithmic problem all orbits 
with nonzero angular momentum are bounded, they do not fall/eject into/from 
the source, and they do not reach maximum distance with respect to the source. 

We return to the analysis of the unperturbed system ()6.2|) and notice that the 
system decouples in the sense that the first two equations are independent of ip. 
Thus, if one solves the equations of r and (p, then ip is obtained by replacing the 
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expressions for r and ip into the third equation and then integrating. Therefore, 
a detailed qualitative analysis of the reduced system: 




(6.8) 



is extremely useful, as it may be extended to the full (r, ip, ip) space by introducing 
the third coordinate ip at the end of our investigation. 

The reduced system is relatively easy to describe. The motion takes place 
on the cylinder [0, R ma x] x S 1 - There are four degenerate saddle equilibria lo- 
cated at (0,0), (0, 7r), (RmaxjO), (Rmax, it), and two centers at (Ro,n/2) and 
(-R ,37r/2), where R is the solution of the equation 2(1 — f(r)) — r 2 = 0. 
A direct computation shows that the eigenvalues for the centers are given by 



Ai j2 = ±ii?oy 2(-Rq + 2)/(Rq + 2). Also, there are four invariant manifolds {r = 
0}, {r = R max }, {(p = 0}, and {(p = n}, forming two heteroclinic cycles connect- 
ing the saddle equilibria (see Figure 3). 

We now lift the dynamics from the (r, ip) phase space into the full (r, ip, ip) 
solid torus by taking into consideration the third ip E S 1 coordinate, as well as the 
dynamics on the collision manifold, zero velocity manifold and the " sin tp = 0" 
zero angular momentum manifolds. 

The global flow, which takes place in the solid torus [0, R max \ x S 1 x §\ can 
be seen in Figure 4 and is represented as follows: 

- the outside boundary corresponds to the collision manifold {r = 0}; 

- the interior boundary corresponds to {r = R max }- The centered circle of the 
torus was "blown up" artificially to an inside torus such that the dynamics on 
the zero velocity manifold can be seen (this is just a visual artifact and does not 
modify the analysis); 

- the motion takes place in between the exterior boundary {r = 0} and the 
interior boundary {r = R ma x}', 

- the surface of zero angular momentum " sin tp = 0" divides the space into two 
distinct global invariant manifolds: one corresponding to ip E (0, 7r) or, equiva- 
lently, to motion with positive angular momentum (C > 0), and one correspond- 
ing to p E (it, 2tt) or simply, to C < 0; the flow is symmetric with respect to the 
horizontal plane C = 0. In this latter plane, the physical motion is rectilinear 
and represents orbits that are ejecting from the collision manifold, reaching the 
maximum distance at the zero velocity manifold and falling back on the collision 
manifold. 

- there are two periodic orbits situated symmetrically with respect to the 
horizontal plane (C = 0), namely: {(r, ip, ip) \ r = Rq, <p = ir/2, ip = ip Q + 
R a s; s > 0, ip — ip(0)} and {(r, ip, ip) \ r = R , ip = 3ir/2, ip = ip Q — R Q s; s > 
0, ipQ = ip(0)}. Around those two periodic orbits, the phase-space is foliated by 
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tori-like surfaces, parametrized by the angular momentum integral ([6.6)1 : 

- the equilibria on the collision manifold are connected with the equilibria on 
the zero velocity manifold, through heteroclinic cycles of the form: 

a — ► b — ► c — ► d — > e — ► / — ► g — ► h — ► a, 

where a = (0,0, ip ), b = (0,ir,if> ), c = (R max ,Tr,ip ), d = (R max ,0,ip + tt), 
e = (0,0,ip + n), f = (0,7r,^ + 7r), 9 = (Rmax, n, ip + vr) and h = (R max ,0,ip + 
2ir) = (Rmax,0,i^o), for each fixed if) Q fixed in [0, 2ir] corresponding one cycle. 

Let us look at the surfaces with positive values of the angular momentum, 
C > (for negative C, by symmetry, the flow is identical but of opposite sense). 
Since each value C = const represents a tori-like surface, the space contains a 
series of invariant manifolds nested one into the other. C varies from to a 
maximal value corresponding to the degenerate torus r = R , i.e. a circle (see 
Figure 5). 

The dynamics will change as one shifts from the tori close to C = towards 
the tori near the periodic orbit r = R Q . In the (ip, if>) plane, the orbits close to the 
collision manifold {r = 0} have very large slopes dip /dip (for example, the orbits 
on cx in Figure 5). The orbits close to the periodic orbit r = R have slopes close 
to —1 (e.g. the orbits on C3 in Figure 5). In between, there is a smooth transition 
in the values of the slopes (one has to keep in mind that the flow is continuous 
and therefore orbits which are initially close have to stay close at all times). 

For a fixed value of C or, in other words, for a fixed torus, the motion is 
bounded between a minimum and a maximum value, which depend directly on 
C. We note that in the real physical space (xi,x 2 ), this translates into a bounded 
motion between a minimum and a maximum value, representing a loop orbit. The 
periodic solutions, (Rq, tt/2, if)) and (Rq, Sti/2, if)), correspond to the two parent 
families of the loop orbits, one counter-clockwise (C > 0) and the other clockwise 
[C < 0) |2] • We shall discuss the family of loop orbits in more detail in Sections 
7 and 8. 

7 Dynamics of the Non Axis- Symmetric System 

In this section, we proceed to investigate the general form of the non axis- 
symmetric system ([6.1)1 . We will limit our analysis to the case in which the 
anisotropy is small (/i ~ 1) and therefore we can treat it as a small perturbation 
(e = JJi — 1 > 0) to the isotropic system. 

In the following, we derive the curves of equilibria on the collision manifold 
and on the zero velocity manifold. We observe that these curves, located in the 
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horizontal plane (C = 0), remain upon the perturbation within 0(e) distance 
of this plane. We then prove that most orbits in the perturbed case avoid the 
origin, a result initially discovered in the numerical study of Miralda-Escude & 
Schwarzschild JT0]- For this, we show that the curves of equilibria on the collision 
manifold are degenerate saddles, i.e. the equilibrium points on these curves 
admit a 2D stable and a 2D unstable manifold. The implication is that, in the 
three-dimensional space, the dimension of the set of initial conditions leading to 
collision is two (i.e. the Lebesgue measure of the set is zero). In physical space 
this translates into a zero probability of finding orbits falling into the source. 

LEMMA 7.1 The equilibria of the vector field (jfi.lj) consists of four closed curves, 
two belonging to the collision manifold and the other two to the zero velocity 
manifold. The collision manifolds curves of equilibria are given by: 

G\ := {(r, (p, if)) | r = 0, ip — arctan[(l — e) tan^>] — if), if) = ifo G [0, 2tt) }U 

U(0, 0,7r/2) U (0,0,37r/2) (7.1) 

and 

C 2 := {(r, (p, if>) | r — 0, (p = tt + arctan[(l — e) tan-0] — tf>, if> = if> G [0, 27r) }U 

U(0,7r,7r/2) U (0,7r,37r/2). (7.2) 
The curves of equilibria on the zero velocity manifold are given by: 

Zi := {(r,ip,if)) | r = R ma x, f = arctan[(l + e)tan^] - if), if) = ip G [0,2tt) }U 

{Rmax, 0, tt/2) U [R max , 0, 3tt/2), (7.3) 

and 

Z 2 := {(r, <p, if)) | r = Rmax-, V 9 = vr + arctan[(l + e) tan^] —if), if) = ipo G [0, 2tt) }U 

(R max , n, tt/2) U (R max , 71, 3vr/2). (7.4) 

PROOF The proof follows by inspection. One can verify directly that the two 
families of curves cancel the vector field. To see that there are no other equilibria, 
suppose r 7^ and r ^ Rmax- Then f = and if) = lead to 

(1 + e) cos(y9 + if)) cos if) = — sin(<yj + if>) sin if) 
(1 + e) cos(y? + if) smif) = sin(<£> + if) cosif). 
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It is easy to check that the above relations cannot coexist. 
OBSERVATION 7.2 d and C 2 are within 0(e) distance of the circles 

{(r,<p,ij>)\r = 0, y? = 0, ij = iJ e[0,2n)} 

and respectively, 

{(r, <p, -0) \r = 0, = 7r, V = V'o e [0, 2vr) }. 

Also, on the zero velocity manifold, Z 1 and Z 2 are within 0(e) distance of the 
circles 

{(r, v?, V>) | r = iCa*, ^ = 0, ^ = e [0, 2tt) } 

and 

{(r, </?, 0) | r = # maa; , (p = tt, ip = ipQ G [0, 2tt) }, 

respectively. 

LEMMA 7.3 The curves of equilibria Z\ and Z 2 on the zero velocity manifold 
are degenerate saddles (see Figure 6). More precisely, for a fixed ipo, the two 
corresponding equilibria 

(R max , arctan[(l + e) tan^o] - ipo, V'o) e ^i 

and 

(-Rmax, 7T + arctan[(l + e) tan^o] - V'o, V'o) £ ^2 
behave like saddles in the (r, </?) plane. 

PROOF Recall from Section 6 that on the zero velocity manifold the flow is 
degenerate and in the (ip, ip) coordinates it reads: 

<P = -Kiaxi^iv + 0) cos ^ - (! + e ) cos (v? + sinijj], 
= Kiaxl^iv + 0) cos^ - (1 + e) cos(v? + ip) simp]. 

By symmetry, it is sufficient to study the flow around one of the curves of 
equilibria, let's say Z\ (around Z 2 we just have to reverse the arrows). For 
simplicity, for a fixed ipo, we denote the ip component of an equilibria by <pz(ipo), 
i.e. (pz(ipo) '■= arctan[(l + e)tan-0 o ] — V'o- Observe that (R m ax,<Pz(i J o),i J o) are 
simple roots for f and (p. Computing the eigenvalues at (R max , ipz(ipo), ipo), we 
obtain: 

A r = ~ or Sma l 2 ( ~ 2f'(Rmax)) [cos p>z(4>o) + e cos(ip z (ipo) + V'o) cosipo] 

A v = -2 cos (ipz(ipo) + V'o) cos^o - esin (tpz(tpo) + ipo) sin^ 
A^ = 0. 
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Recall that if) is fixed in [0, 7r/2]. If ipo = tt/2, we obtain (pz (tpo) = < / 7 ( 7r /2) = 
and 



p3 

^ mttc j 



-2f'(R max ) > 0, 



A v = -e < 0, 
^ = 0. 

Therefore (R ma ,x, 0, 7r/2) is a degenerate saddle. If V'o £ [0, 7r/2), we have: 



A r = -(1 + e)- 



^Li + 2 
A^, = -(2 cos 2 *0o + e sin 2 V'o) 



2f(Rmax)) 



COS 



(y?zOo) + ^o) 



COS l[)Q 

cos ((pz(i>o) + i/>o) 



COS '0, 







Since tp G [0, 7r/2) and (pz(i>o) = arctan[(l + e)tan-^o] — i>o, it results that 
fz{^o) +^o£ [0, 7r/2) and furthermore, cos {ipz(tpo) + ^>o) > 0. Therefore A r > 
and A^ < 0. 

The other cases where ip G (■k/2,tt) U [7r, 37r/2] U (37r/2,27r) can be treated 
similarly, reaching the same conclusion, i.e. the equilibria points Z 1 are all de- 
generate saddles, with: 

' \ r > 0, 

\ < o, 

A^ = 0. 



LEMMA 7.4 The curves of equilibria C\ and C 2 on the collision manifold are 
degenerate saddles (see Figure 6). More precisely, for a fixed ip , the two corre- 
sponding equilibria 

(0, arctan[(l — e) tan V>o] _ ipo, ipo), 

(0, 7r + arctan[(l — e) tan^o] _ i>o, i>o) 
behave like degenerate saddles. 

PROOF The equations of the flow on the collision manifold are 

Lf> = 2[(1 + e) sin((/9 + ip) cosip — cos(</? + ■?/>) sin^] 
ip = 

As in the previous proof, it is sufficient to investigate the flow around the 
equilibria curve C\ (around C2, by symmetry, we just have to reverse the ar- 
rows). Fixing ipQ G [0, 2ir), we denote the corresponding fixed point on C\ by 
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(0, (pc(if>a), ipo) and proceed to calculate its eigenvalues. We obtain: 



A r = 0, 
A^ = 2 
A^ = 0. 



(1 + e) cos (y?cOo) + ipo) cosipo + sin (<pc(ipo) + ^o) sin^ 



We observe that the vector field ()6.1|) manifests a degeneracy around r = 0. This 
means that any fixed point on the curve C\ might have a non-hyperbolic character 
in the (r, if) plane. For the moment, let us discuss the sign of A^,. Similar to the 
analysis in the case of the curve Zi, we fix ip £ [0,7r/2]. If ipo = vr/2, we have 
Pci 11 /2) = and therefore A^ = 2 > 0. If ip G [0, vr/2), we obtain: 

sin(^o) 
\ 



(<Pc(ipo) + ipo) 

Since (<pc(ipo)+ipo) = arctan [(1— e) tan^o] G (0, vr/2), it follows that sin((^c*('0o) + 
ipo) > and therefore A v > 0. The same type of reasoning applies for all the other 
cases, i.e. for ip Q G (vr/2,7r) U [vr,37r/2] U (3vr/2,27r). Therefore we have proved 
that A^ > for any fixed ipo G [0, 27r). 

In conclusion, the flow around the equilibria (0, (pc(ipo), ipo) G Ci is degener- 
ate, with eigenvalues: 

'x r = 
\ > o 

A^ = 0. 

Obviously, the linear approximation of the flow does not provide enough in- 
formation about the behavior around the fixed points. Let us take a closer look 
at the vector field around (0, (pc{ipo), ipo) £ Ci- 

Notice that on the ip direction the flow is null, as every point on the circle 
ip G [0, 2vr) is a parameter for a fixed point. It remains that in order to describe 
the asymptotical behavior around (0, <pc(ipo), ipo), one has to investigate the flow 
in the (r, (p) coordinates. A, r = generates the center manifold E c -the span 
of the zero eigenvector ([7], PH]). The general theory ensures the existence of 
an invariant manifold W c tangent to E c at (r, (fc{ipo))- W c may not be unique 
and, usually, it involves a loss of smoothness. Also, around the fixed point, W c is 
described by a one-parameter family of curves, i.e. W c = {(r, p) \ ip = pp(r), [3 G 
E}. 

In our case, we will compute W c directly near (r, pc{ipo))i as follows (see 
for example, ;.8]): We know that by a proper transformation of coordinates the 
equations for (r, p) have the structure 



j. 3 



r 



-2(1 -/(r))F!(^V;e 



+ 2 v M /y nr,r ' y (7.5) 



ip = X^ip - (p c (ipo)) + 
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Near r = 0, we can make the approximation: 

r 3 „, , s ,1 



2 2 2(1 - /(r)) = r 3 -(l + r'/2 + . . .)2(1 - /(r)) = (7.6) 

= r 3 + higher order terms 

(note that we are not expanding /(r), a function that is only differentiable, but 
we are merely looking for the dominant term as r goes to zero). Therefore, around 

(r, <fc(i>o)) we have: 

r = -ar 3 + . . . , 
= - <Pc(il>o)) + ■ ■ ■ , 
where a := Fi(<pc(ipo),ipo', e) is a positive number. We obtain immediately that: 

dr (7.8) 



with solutions: 

, v / Pe x - /{2r2) if r > q\ 
^ (r) = 1 if r = 0, (7 ' 9) 

where f3 G R is a parameter. 

We sketch the family of solutions ^pp{r) in Figure 7. As it can be easily seen, 
the fixed point (0,(pc(ipo)) is indeed a saddle. 



Since the above reasoning applies for any of the points of equilibria, and using 
the Lemmas 7.1 and 7.3, we can state the following: 

THEOREM 7.5 Let be the system (|6.1|) . Then the equilibria are given by the 
curves Ci, C2 Z\ and Z2, as defined in Lemma 7.1, and each of this curves admits 
a two-dimensional stable manifold and a two-dimensional unstable manifold. 

COROLLARY 7.6 In the anisotropic logarithmic problem the Lebesgue measure 
of the set of initial conditions that lead to collision is zero. 

The last result states that in the phase space the probability of choosing any 
initial conditions that lead to collision is zero. The set of such initial conditions 
is formed by a two-dimensional manifold embedded in a three-dimensional phase 
space. Since we reach a similar conclusion in the isotropic case, we can conclude 
that for the singular logarithmic potential the anisotropy does not increase the 
probability of finding orbits falling into the source. In other words, we have proved 
analytically that in the non axis-symmetric case all orbits (with the exception 
of a set with measure zero) are centrophobic, a result that was also noticed in 
previous numerical studies [TP] . 
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We turn our attention now to the behavior of the angular momentum C(s). 
Recall that C(s) was given in (|4.9|) . For h = and in terms of e and (r,ip,if)), 
the relation (14.911 becomes: 



C(s) = g(r)2[l — f(r)] (1 + e) sin (ip + ifA cos if) — cos(y? + ip\ sin if) . (7.10) 

We note that under the perturbation, the horizontal surface C = (see Figure 
4) splits up, and some of its vestiges are to be found along the two-dimensional 
stable manifold of C% and along the unstable manifold of Ci- In physical space, 
this corresponds to the case of rectilinear orbits (C=0), which are falling into or 
ejecting from the source. 

The variation of angular momentum 1)4.14)1 reads in our coordinates: 

dC 

— = -eg(r)2[l - f (r)} 3/ 2 simp cos ip. (7.11) 

Since the derivative of C is bounded, there are no "blow-up" type effects in 
the evolution of C. Also, since the product g(r) ■ 2[1 — f(r)) is always positive, it 
follows that the critical points of C(s) correspond to if) G {0, tc/2, it, 37r/2} (see 
Figure 8). These orbits, for which C(s) displays a sinusoidal- type behavior and 
admits four critical points, represent the family of loop orbits (see also a similar 
result in the study of Touma & Tremaine ) . 

On the other hand, there are orbits which, under perturbation, will slip in 
between the stable manifold of C\ and the unstable manifold of C2, switching 
the sign of the angular momentum. By the nature of the phase-space, which is 
a solid torus, these orbits must wind around indefinitely. These orbits pertain to 
the family of "boxlet" orbits discovered by Miralda-Escude & Schwarzschild |TU] . 
For large perturbations, more and more orbits will break away from the curves of 
equilibria C\ and C2, and become boxlet, a result which was also pointed out in 
[TO] • The separatrix that divides the phase-space between loop and boxlet orbits 
intersects the (C = 0) plane at if) = n/2 and if) = 3n/2 (see [H3). Some of the 
orbits that wind up inside the torus will eventually close, becoming the parent 
orbits of the resonant families (i.e. banana, fish, pretzel, etc) 



We do not present here a rigorous proof for the existence of the boxlet obits or 
the resonances, leaving it for a future study. However, we present below a partial 
result concerning the orbits which preserve the sign of the angular momentum. 

Let us notice that the angular momentum can be regarded as a function of 
two arguments, namely as C = C(s, e). Then, around the equilibrium point (s, 0), 
we have the following approximation: 



C(s,e) = C(s,0) + e dC{S > e) 



+ 0(e 2 ) (7.12) 
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or, in our coordinates, 

C(s, e) = C(s, 0) + eg{r)2[l - /(r)] sin (<p + ifj) cosip. (7.13) 

But C(s, 0) is a constant, since this is the case of the unperturbed motion, where 
the angular momentum is an integral of motion. Denoting Co := C(s,0), it 
follows that: 

C(s, e) = C + eg{r)2[l - /(r)] sin(y? + if) cosip. (7.14) 

Since eg(r)2[l — /(r)] sin (yj + ?/>) cos^ is bounded, we conclude that C(s) pre- 
serves its sign for large values of Co and small e. This is best noticed near the 
unperturbed periodic orbits, for example, (r,<p,ip) = (Ro, ir/2, ip), ip G [0,27r), 
where Co = e _1 ^ R °2(l — R^\yiRq) and where the perturbation induces loop or- 
bits (see Figure 9). 



8 The Orbital Structure 

The orbital structure can also be retrieved numerically, by integrating the general 
form of the system (|6.1|) . The return to the initial coordinates in the physical 
space, (xi,x 2 ), can then be made through the relations: 

xi = qi/fi= • re ( ~ 1/r2) cos(y? + ip) 

%2 — Q2 — re^~ x l r ' sin(</? + ip), 

and recalling that the time scale was modified such that both the singularity and 
the zero velocity manifold are now reached in an infinite time. 

By appropriately choosing the initial conditions and recalling that e = yjji — 
1 = 1/6 — 1, one can retrieve the resonance families. For example, we know that 
the family of loop orbits (1:1) develops around the two periodic orbits r = R 
(see Figure 4), i.e. near: 

x x = (l//i) • re^ 1 ^ cos(tt/2 + ip), 

x 2 =re(- 1 / r2 )sin(7r/2 + V), (8.2) 
ip = R%s + ip . 

and 

x x = (l//i) • re^" 1 ^ 2 ) cos(3tt/2 + ip), 

X2 = re (-i/r 2 ) s in(3vr/2 + 0), (8.3) 
ip = -R 2 s + ip . 

The two periodic orbits correspond to the two parent orbits of the loop family, 
one evolving clockwise (C < 0), the other counter-clockwise (C > 0). Integrating 
the system (|6.1|) for the case e = 0, around either (Ro,ir/2,ip) or (Ro,37i/2,ip), 
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one retrieves a circle of radius Ro, as indicated with thick line in Figure 9. The 
same figure shows, with dotted line, the orbital structure in the vicinity of the 
parent orbit for the case when e = 0.1 (corresponding to b = 0.9). 

Besides the loop resonance, several other resonances can be recovered from 
the new system of differential equations. We remark that under large pertur- 
bations (i.e. large deviations from axis-symmetry), the families of minor orbits 
will occupy more and more of the phase-space. We do not attempt here to cover 
the entire phase-space, but rather give a few examples of families of minor orbits 
found with our new system of equations. For example, Figure 10 shows an exam- 
ple of an orbit around the fish (3:2) resonance, in the case e = 0.3 (corresponding 
to b = 0.7) and the initial conditions r = 0.9, <p = n/12 and if) = tt/8. Figure 11 
shows an orbit around the pretzel (4:3) resonance, when e = 0.1 and r = 0.45, 
(p — 0, if) — 7r/ 14. Figure 12 shows an orbit near the (5:3) resonance, obtained for 
e = 0.3, r = 0.4, ip = 7r/12, if) = 7r/6. We note however, that these parameters 
should be taken only as a guidance for the location of the resonances. 

9 Conclusions 

Previous numerical studies have revealed several important aspects of the orbital 
structure of the logarithmic potential: the division between loop and box orbits, 
the presence of resonances, the scattering effect of the singularity (which renders 
the box orbits unstable), and the transition to chaos. We have performed an 
analytical study of the singular logarithmic potential and proved several of these 
results. 

We summarize our results as follows: 

- we provide a description of the dynamics near the singularity and at the 
maximum distance from the source permitted for a given level of energy. 

- in the axis-symmetric case, we retrieve the complete global dynamics of the 
orbits and describe it on a solid torus bounded by the two surfaces {r = 0} 
and {r = R max }. We find analytically the two periodic orbits r = R , which 
correspond in physical space to the two parent families of the loop orbits (in 
clockwise and respectively, in counter-clockwise motion). 

- in the non axis-symmetric case, we prove that all orbits, except a negligible 
set, are centrophobic -a result that has been originally discovered in the numerical 
study of Miralda-Escude & Schwarzschild [TU] . 

- in the same non axis-symmetric case, we show that there exist orbits which 
preserve the sign of the angular momentum and retrieve the loop resonance. 
Finally, we also show how several other minor family orbits can be obtained from 
the our new system of equations. 

The analytical description in general non axis-symmetric case remains still 
open, several problems requiring further investigation. One of them is the re- 
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trieval of the family of boxlet orbits, which are known to dominate the dynamics 
near the singularity [TU]. We conclude, however, that under large perturbations 
(large deviations from axis-symmetry), most of the phase space is occupied by 
orbits that slip in between the stable manifold of C\ and the unstable manifold 
of C2 and wind around indefinitely. Most of these orbits pertain to the family of 
the boxlet orbits. Some of these winding orbits may be closed and would lie at 
the origin of the resonant families observed experimentally. Further work in this 
direction will have to include an in-depth analysis of resonances in the (21,0:2) 
variables and a Fourier expansion of the periodic solutions for r(s). 

Another aspect that remains to be clarified is the existence of the stochastic 
orbits near the singularity and the transition to chaos. The analytical study in 
this case is hindered by the fact that the system given by the Hamiltonian (j2.1j) 
does not admit any hyperbolic equilibrium points, and therefore, the perturbative 
methods which are usually employed in proving the existence of chaos (including 
the Melnikov method [H]) become unapplicable. We note that the Melnikov 
method has been used in the past in the case of the logarithmic potential 5j, 
however not on the exact Hamiltonian, but on the integrable Stackel Hamiltonian 
|12j . Understanding the onset of chaos in the singular logarithmic potential is a 
very important problem for the construction of galaxy models based on libraries 
of orbits (see for example, |14j ) - in which generally, it is a priori assumed that 
the stochastic orbits play a negligible role. 
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Figure 1: On the collision manifold, the flow follows the parallel lines ip = V'o = 
const. The flow vanishes on the equilibrium curves C + and C~. 

Figure 2: On the zero velocity manifold, the flow follows the parallel lines 6 = 6q = 
const. The flow vanishes on the equilibrium curves V + and V~ . 

Figure 3: The reduced phase-space (r, if) in the case fi = 1. There are six equilibrium 
points: two centers, (Rq,it/2) and (i?o,37r/2), and four saddle equilibria, (0,0), (0,7r), 

(Rmax, 0), (R max > " ) ■ 

Figure 4: The global flow in the axis-symmetric case (p, = 1), described in a solid 
torus delineated by the two surfaces {r = 0} and {r = R ma x}- The full lines denote 
the heteroclinic cycle for the case Vo = 0. The horizontal plane (C = 0) divides the 
phase-space in two symmetric invariant subspaces (see text for details). 

Figure 5: Section on the Solid Torus at ip = ipQ. Around the two periodic orbits 
r = Rq, the phase-space is foliated by tori- like surfaces of constant angular momentum, 
C. 

Figure 6: The curves of equilibria Z\ and C\, located on the collision manifold and 
on the zero-velocity manifold, respectively. Each equilibrium point on these two curves 
behaves like a degenerate saddle. 

Figure 7: Curves in the y>p(r) family around the saddle point (0, ipc(ipo)). 

Figure 8: The curves of equilibria Ci, C2, Z±, Z2 under the perturbation e. The 
arrows show the behavior of the flow in their neighborhood. The manifolds C\ and C2 
intersect transversely along the lines sintp = and cos ip = 0. 

Figure 9: The circle r = Rq (thick line) represents one of the parent families of 
loop orbits. The dotted line shows the orbital structure that is triggered by applying 
a perturbation (e = 0.1), to the periodic orbit (shown here is the integration only over 
a finite time). 

Figure 10: Orbit around the fish (3:2) resonance. 

Figure 11: Orbit around the pretzel (4:3) resonance. 

Figure 12: Orbit around the (5:3) resonance. 
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